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ABSTRACT 

The validity of the hypothesis that the massive black holes in high redshift quasars grew from stellar-sized 
"seeds" is contingent on a seed's ability to double its mass every few ten million years. This requires that the 
seed accrete at approximately the Eddington-limited rate. In the specific case of radiatively efficient quasiradial 
accretion in a metal-poor protogalactic medium, for which the Bondi accretion rate is often prescribed in 
cosmological simulations of massive black hole formation, we examine the effects of the radiation emitted 
near the black hole's event horizon on the structure of the surrounding gas flow. We find that photoheating 
and radiation pressure from photoionization significantly reduce the steady-state accretion rate and potentially 
render the quasiradial accretion flow unsteady and inefficient. The time-averaged accretion rate is always 
a small fraction of the "Bondi" accretion rate calculated ignoring radiative feedback. The pressure of Lya 
photons trapped near the H II region surrounding the black hole may further attenuate the inflow. These results 
suggest that an alternative to quasiradial, radiatively efficient Bondi-like accretion should be sought to explain 
the rapid growth of quasar-progenitor seed black holes. 

Subject headings: black hole physics — cosmology: theory — galaxies: active — galaxies: formation — 
hydrodynamics — quasars: general — radiation mechanisms 



1. INTRODUCTION 

The origin and the early growth of massive black holes re- 
mains poorly understood. The massive black holes in quasars 
(M BH ~ 10 8 - 1O 1O M ), active galactic nuclei (AGN; M BH > 
1O 5 M ), and the quiescent nuclei of nearby galaxies may 
have started out as stellar-mass (Mbh < 100 Mq) "seed" black 
holes (e.g., |Madau & Rees||2001| |Menou etaT1|2001| |Islam 
2003| >. Is this plausible, that is, could the seed black 
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holes have grown rapidly enough in the cosmic time avail- 
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nd references therein)? The 
rate at which a seed massive black hole can accrete is lim- 
ited by the local density and the thermal structure of the pro- 
togalactic medium and by the effects of the radiation emit- 
ted near the event horizon on the accretion flow. Cosmo- 
logical hydrodynamic simulations suggest that gravitational 
collapse produces dense central gas concentrations in proto- 
galaxies (e. g., |Springel et al.|2005||Li et al.|2007||Pelupessy| 
et al.|20071 Wise & Abel 200 7a|bHWise et al.|2008HDi Mat-| 
teo et al. 2008 1 Greif et al. 2008 ). Atomic densities have been 



10 4 



cm" J (e.g., |Greif et al.|2008| , and as 



found to reach n ■ 

much as n ~ 10 6 cm" 3 averaged over the central parsec around 



the potential minirrmm ( Bromm & Loeb|2003| |Wise & Abel 
|2007b[|Wise et al.|2008 1. On spatial scales that are resolved 
in the simulations, gas is sufficiently concentrated to enable 
rapid accretion onto a seed black hole. An exception are the 
first hundred million years after the seed was formed, during 
which the surrounding gas density is lowered by the radia- 
tive feedback fr om the bl ack hole's progenito r star (see, e.g., 
| Johnson & Bromm|2007l |Alvarez et al.|2006||2008] l. 

Given an ample gas supply, will rapid accretion be inhibited 
by radiative effects? A reassessment of an accreting black 
hole's ability to control its own gas supply is needed to im- 
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prove the realism of the treatment of black hole accretion in 
cosmological simulations. Existing cosmological simulations 
modeling the growth of seed black holes do not resolve the 
spatial scales on which some of the radiative processes may 
alter the accretion. The simulations also do not resolve the 
fine structure of the dense, turbulent, and possibly multiphase 
protogalactic medium in which the black holes are embed- 
ded. Semianalytic prescriptions are normally adopted for the 
accretion rate, but these prescriptions normally do not take 
into account the radiative feedback; it is normally assumed 
that a black hole, residing in a pressure-supported primordial 
gas cloud, can accrete steadily at the Bondi rate subject to 
the Eddington limit (e.g., Volonteri & Rees|2005 Alvarez et 
al.|2006||Johnson & Bromm|2007[ |Pelupessy et al.|2007| |Di 



Matteo et al.|2008[|Greif et al.|2008| l. Here, we will evaluate 
the applicability of this assumption in view of the local radia- 
tive feedback that is present if the black hole accretes in a ra- 
diatively efficient fashion. We restrict our analysis to the early 
growth of quasar-progenitor "seed" black holes, which are oc- 
casionally referred to as "miniquasars," where a stellar mass 
or an intermediate-mass black hole (1O 2 M <M BH < 1O 5 M ) 
accretes from a metal and dust poor environment. 

The black hole's growth rate is particularly sensitive to the 
detailed thermal state of the irradiated accretion flow. This 
can be seen by noticing that the accretion rate, for quasiradial 
accretion, is influenced by the conditions at the sonic radius 



r s ~3x 10 — cm, 

^s,5 



(1) 



where r s = 10 5 r s 5 K is the temperature of the photoion- 
ized and photoheated flow at the sonic radius and Mbh = 
100 MjMq is the black hole mass. The sonic radius is nor- 
mally unresolved in cosmological simulations of accretion 
onto black holes in the intermediate range of masses. 1 It was 

1 The recent, highly- resolved simulation of primordial protostar formation 
by |Yoshi da et al. 1 2008 1, which has the requisite spatial resolution to resolve a 
some radius it it exists, only proceeds to the point where the initial hydrostatic 
core is formed, and does not treat the subsequent accretion flow onto the 
growing core. 
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recognized early that photoheating and photoionization pres- 
sure may prohibit steady radiatively efficient accretion (e.g., 
Shvartsman|1971||Buff & McCray|1974j|Hatchett et al.|1976[ 
Ostriker et al. 1976). It was suggested that accretion can still 
proceed at rates disallowed by the steady state solutions if 
cycles of rapid gas inflow, during pauses in accretion near 
the event horizon, alternate with photoheating or photoioniza- 
tion pressure-driven outflo ws ( Buff & M cCray 1 1 974| |Ostriker| 
et al.||1976| |Cowie et al.|[l978| |Stellingwerf & Buff|[T982f 
Begelman 1985). Such quasiperiodic cycling is seen in one- 
dimensional simulations of Compton-heated accretion onto 
Mbh <; 1O 8 M black holes in galaxy clusters, where the black 
hole accretes from a hot, ionized, and pressure supported at 



UV and X-ray radiation is sensitive to metal abundances, es- 
pecially for T < 10 4 K (e.g.,|D onahue & Sh ull|l9 9lt and for 
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mosp here ( |Ciotti & Ostriker|19 97 2001 2007] |Sazonov et al. 
|2005) l. Recently, |Ricotti et al.| ( |2008| > revisited the problem 
of irradiated quasiradial accretion in the context of primordial 
black hole growth following the cosmic recombination (see, 
e.g., Ricotti 2007, and references therein), and suggested that 
the accretion duty cycle is determined by the periodic forma- 
tion of an H II region surrounding the black hole. 

It was further recognized that the formal existence of 
steady-state, spherically symmetric accretion solutions is sen- 
sitive to the treatment of boundary conditions far from the 
sonic radius (e.g., Bisnovatyi-Kogan & Blinnikov|1980" i, and 
that these accretion flows can be locally thermally unstable 
(e.g., |Stelhngwerf|1982[|Krolik & London|1983| l and should 
break down into time-dependent two-phase structure, contain- 
ing a warm ionized phase and a hot, coronal phase (e.g.JKro- 



lik et al.|[T9"8"T] >. |Wang et al\ ( |2006| ) claimed that Compton 
heating in the vicinity of a seed massive black hole reduces 
the radial accretion rate to a small fraction of the Eddington- 
limited rate; we here suggest, however, that thermal runaway 
may engender the Compton-heated coronal phase only in 
metal-rich flows, where photoionization heating of the incom- 
pletely stripped oxygen drives gas heating beyond ~ 10 5 K 
(see, e.g., |Kallman & McCray fT982| and §[XT]below). 

The impact of the radiation field produced near the event 
horizon on the accretion flow and on the state of the in- 
terstellar medium of the protogalaxy and that of the inter- 
galactic medium (e.g., Dijkstra et al. 2004; Kuhlen & Madau 
2005) |Zaroubi et al.||2007| |Thomas & Zaroubi||2008| |Ripa-| 
monti et al.||20081 ISpaans & Meijerink||2008jiis sensitive to 



the shape of the spectral energy distribution (SED) of the cen- 
tral source. The SEDs of rapidly accreting low-mass massive 
black holes (M B h ~ 10 s - 1O 6 M ) exhibit significantly larger 
X-ray (2 keV) to optical spectral ratios than the AGN con- 
taining more massive rapidly accreting black holes ( Greene & 
|Ho|2007] l, as is expected if a fraction of the radiation is pro- 
duced in a geometrically thin disk. On the low mass end, if the 
microquasar SEDs (e.g., Remillard & McClintock 2006] and 
references therein) are an adequate prototype, the seed black 
hole SEDs may contain energetically significant components 
extending into the hard X-rays. The luminosity-weighted av- 
erage spectrum of AGN containing intermediate mass black 
holes could differ substantially (see, e.g., Venkatesan et al. 



2001 1 |Madau et al.||2004[) from the scaled average quasar 



spectrum of Sa zonov et al.| p004). The ability of X-rays to 
escape the protogalaxy affects their contribution to the soft 
X-ray background (e.g., Venkates an et al.|2001| Dijkstra et 
al.|2004 Salvaterra et al.|2005| > and the infrared background 
( |Cooray & Yosh ida 2004). Another difference between the 
first AGN and the starburst or post-starburst AGN is related 
to the differences in metallicity of the accretion flows. The 
thermal phase structure of the interstellar medium exposed to 



T > few x 10 4 K (e.g., |Kallman & McCray|1982| and §_ 

The role of radiative feedback in the formation of the first 
massive black holes resembles the radiative regulation in the 



formation of the first massive protostars ((Omukai & Palla 


2001 


2003; Omukai & Inutsuka 2002; Omukai & Yoshii 


2003 


Tan & McKee 2004; 


McKee & Tan 2008 ), though, of 



course, the central sources have very different spectra. In pro- 
tostars, an H II region forms around the protostar and the pro- 
tostellar disk that feeds its growth; persistent accretion onto 
the protostar may be quenched by radiation pressure. The 
protostellar accretion is characterized by a lower radiative ef- 
ficiency and a higher accretion rate than the black hole ac- 
cretion for the same accretor mass. The growing protostar is 
embedded in a supersonically collapsing and very dense en- 
velope from inception, whereas here we assume that the seed 
black hole is born in the collapse of a massive star without 
such an envelope. The presence of an infalling envelope im- 
plies that the H II region is initially trapped near the protostar, 
inside the radius where the envelope infall velocity turns from 
subsonic to supersonic. For accretion from a diffuse medium 
onto a seed black hole, the H II region is inevitably extended 
and the flow crossing the ionization front is highly subsonic; 
this marks a crucial difference with the protostellar accretion 
scenarios. 

Ignoring radiative effects, accretion onto the black hole is 
quasiradial on certain length scales (e.g., r s < r < 10 pc) if tur- 
bulence in the gas is weak (Krum holz et al.||2006| l and if the 
gas is not rotationally supported on these scales. If the accre- 
tion flow is shock-free (an unlikely condition) and possesses 
small net rotation, the buildup of vorticity near r ~ r s may 
reduce the accretion rate by ~ 60% (Krum holz et al.| ;2005). 
Quasiradial accretion may be expected even when the baryons 
in a protogalaxy initially form a rotating disk, because self- 
gravity in the gas on scales of the protogalactic disk destabi- 
lizes rotational equilibria to convert disk-like configurations 
into quasiradial, stratified, pressure supported, and possibly 
turbulent configurations. The gas distribution could still be 
rotationally supported on larger scales, where, e.g., dark mat- 
ter dominates gravity, and on much smaller scales, where the 
black hole dominates gravity. This description may apply to 
high-redshift protogalaxies (see, e.g., Oh & Haiman|[2002{ 



lyblonteri & Rees 2005; Wan g et al.||2006| ), and scThere, we 
focus on angular momentum-free accretion and defer exam- 
ining the role of angular momentum to a subsequent paper. 

The quasiradial gas flow can either be steady and radial, or 
unsteady and characterized by alternating inflow, outflow, and 
nonradial motions. In view of these possibilities, this work is 
organized as follows. In §[2] we attempt, and fail, to construct 
a steady, radial solution for accretion at high accretion rates 
and high radiative efficiencies. In § 3]we provide a qualitative 
analysis of time-dependent, episodic accretion, and attempt 
to estimate the average accretion rate. In §|4]we discuss the 
consequences of the presence of cold, inhomogeneous, and 
turbulent gas in the vicinity of the black hole. In §|5] summa- 
rize our main conclusions. In Table [T] we present an overview 
of our notation. 

2. THE PROSPECT OF TIME-INDEPENDENT ACCRETION 

In this section we set out to test the model, ubiquitous in 
semi-analytic and semi-numerical studies of massive black 
hole evolution, in which the accretion onto a black hole from 
the protogalactic medium is steady and quasiradial. In § 2.1 
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Quantity 

Acceleration due to gravity 
Acceleration due to radiation pressure 
Spectral index of the SED 
Case B recombination rate for hydrogen 
Collisional ionization rate for hydrogen 
Abundance of species i relative to hydrogen 
Logarithmic slope of r s as a function of e 
Radiative efficiency 

Critical efficiency for radiation pressure suppression at low efficiencies 
Ionization potential of species i 
Duty cycle for episodic accretion 

Density enhancement in the ionized gas in episodic accretion 
Fraction of photon energy going to photoionizations in the ionized gas 
Fraction of photon energy reprocessed to hya 
Fraction of L that reaches the edge of the H II region 
Pressure enhancement due to turbulence in the neutral gas 
Adiabatic index 

Luminosity in units of the Eddington luminosity 

Luminosity of isothermal accretion ignoring radiative heating 

Lya-pressure-limited accretion rate 

Peak luminosity in episodic accretion 

Luminosity of steady-state photoheated accretion 

Luminosity emitted at r < 

Eddington luminosity for Thomson scattering 

Mean molecular mass of the ionized gas 

Turbulent Mach number of the neutral gas 

Mass of the black hole 

Central accretion rate 

Central accretion rate ignoring radiation 

Density at the H II region's edge in a protogalactic density cusp 

Density of the neutral gas and ambient density 

Particle density within the H II region 

Maximum density for high ionization at the sonic radius 

FT-dissociating photon production rate 

Lyman-Werner photon production rate 

Total ionization rate in the H II region 

Total recombination rate in the H II region 

Number of Lya reflections on H II region's walls to escape 

Dimensionless photoionization pressure acceleration in the H II region 

Dimensionless radiation pressure acceleration at the H II region's edge 

Pressure of the neutral gas 

Pressure of the Lya radiation 

Radius of the H II region in a protogalactic density cusp 
Bondi radius ignoring radiative heating and acceleration 
Sonic radius 
Disk radius 

Radius where thermal time equals inflow time 

Radius of the H II region 

Width of the Fb photodissociation shell 

Width of the neutral shell surrounding the H II region 

Line-center Lya scattering cross section 

Flux-averaged photoionization cross section of species r' 

Photon number-averaged photoionization cross section 

H2 formation suppression factor due to FT photodissociation 

Lya line-center optical depth of the neutral shell 

Bremsstrahlung cooling time 

Compton heating time 

Cooling time 

Heating time 

Photoionization heating time 
Inflow time at the sonic radius 
Salpeter mass-exponentiation time scale 
Temperature of the ambient neutral gas 
Temperature within the H II region 
Temperature at the sonic radius 

Density enhancement over isothermal accretion with 7hh = T 5 

Radial inflow velocity 

Dimensionless inflow velocity 

Ionization parameter 

Dimensionless ionization parameter 

Dimensionless radius 

Metallicity 
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we review the standard theory of photoionized radial accre- 
tion, and pay particular attention to the dependence of the 
structure of the flow on the metallicity of the accreting gas. 
We find that at low metallicities, the flow evades thermal run- 
away and heating to the Com pton temperature at least until it 
passes the sonic radius. In § |2.2| we estimate the size of the 
H II region surrounding the black hole and derive conditions 
under which the self-gravity within the ionized sphere can be 
ignored. In § |2.3| we take a closer look at the state of the gas 
as it passes the sonic radius and check whether it is in local 
thermal and statistical equilibrium. In § 2.4 we justify our ref- 
erence choice for the ambient density in the neutral medium 
surrounding the H II region. This justification is necessary be- 
cause the central gravitational collapse in a protogalaxy will 
yield a wide range of densities, yet most of our estimates de- 
pend on a specific choice of density. In § |2.5| we study the ef- 
fects of the photoionization radiation pressure within the H II 
region on the structure of the steady-state accretion flow. We 
find that, depending on the parameters, photoionization radi- 
ation pressure in the outer parts of the H II region may pre- 
vent accretion at near the Eddington-limited rate. In § 2.6 we 



briefly address the case of radiatively-inefficient accretion. In 
§ |2.7| we estimate the pressure of Lya line radiation produced 
in and near the H II region and confined by resonance line 
scattering in its vicinity. We find that Lya radiation pressure 
can exceed the thermal gas pressure in the H II region, and this 
presents an additional challenge to strictly stationary quasira- 
dial solutions. 

2. 1 . Photoionized Quasiradial Accretion 

We ignore the angular momentum of the gas, which is as- 
sumed to be free of metals and dust, on radial length scales 
> r s and assume that the accretion flow becomes rotationally 
supported and collapses into a hypothetical geometrically thin 
disk only at radii < r^k <C r s . Furthermore, we assume that 
the disk accretes onto the black hole with a high radiative ef- 
ficiency e, as is expected for thin-disk accretion, such that the 
bolometric outward radiation flux passing through radius r s is 
F(fs) = eM(r s )c 2 /47rr 2 . We assume that, absent feedback ef- 
fects, the density scale-height is much larger than any other 
length scale under consideration, so that the ambient medium 
has effectively constant density. We restrict our attention to 
the accretion flow at radii r > r s , where we assume that the 
flow is quasiradial and exposed to the radiation emitted by the 
disk at r ~ 0. Without angular momentum, the radial support 
against the black hole's gravity must arise from gas pressure 
gradients and the radiation pressure force. Conservation laws 



yield the relation (see, e.g., |Ostriker et al. 
|Cassinelli|1999) 



1976 Lamers & 



dv 
dr 



jkT 
fi m p 



■■y\ ^^+atot )+(7-l)(#-C). (2) 

r fi nip 



where v is the radial inflow velocity, which is positive when 
the flow is directed inward, 7 the adiabatic index that defines 
the relation of internal energy density to pressure, k the Boltz- 
mann constant, [i the mean molecular mass, m p the proton 
mass, T the gas temperature, a tot = a ra d + a gra v the sum of the 
accelerations due to gravity and radiation pressure, H the pho- 
toheating rate, and C the cooling rate (both per unit mass). 

there exists a radius r, 

W/(7 



As we show in § 2.2 
which the heating time fh ea t 



eq m outside of 
l)fim p H and cooling 



gas are generally short, outside this radius the gas reaches an 
approximate local thermal and statistical equilibrium, and so 
the last term in equation Q that is proportional to H — C can 



be dropped. An estimate that we provide at the end of § 2.2 



below suggests that for metal-poor accretion, the radius r equ ; 
is usually at most slightly larger than the radius r s at which 
the flow becomes supersonic. We ignore this complication 
and assume r equ ; < r s . 

The flow is not adiabatic, but we let c s = (7W / ' [im p ) x l 2 de- 
note the usual adiabatic sound speed. Continuity at the sonic 
radius r s where v = c s requires that 2[c s (r s )] 2 r~ 1 +a grav (r s ) + 
<Zrad( r s) = 0, and because the self gravity of the gas is neg- 
ligible, Ggrav^s) = -GMbh/^. Since the gas is almost fully 
ionized at r ~ r s , the radiation pressure at the sonic ra- 
dius is mainly due to electron scattering. Then a ra d(r s ) = 
-£(r s )a grav (r s ), where £{r) = L(r)/LEdd is the ratio of the to- 
tal luminosity to the Eddington luminosity for all opacities. 
Since under a wide range of conditions (see, e.g., |Blondin| 
1986), the flow at the sonic radius is optically thin to elec- 



time f coo i ~ kT/Qy— \ )nm p C are much shorter than the inflow 
time ~ r/v. Since photochemical time scales in dense, ionized 



tron scattering, provided that the small accretion disk inside 
the sonic radius does not shadow and reprocess to low fre- 
quencies a substantial fraction of the central luminosity, we 
can assume that L(r s ) = L(Q). Then, in a steady state, the lu- 
minosity can be related to the total mass flux into the central 
source, L(0) = eMc 2 = 47rer 2 c s (r s )/i«(r s )m / ,c 2 , where n(r) is 
the gas number density. 

Let H = L/4irr 2 nkTc denote the dimensionless ionization 
parameter introduced in Krolik et al. ( 1981} , which is the ra- 
tio of the radiative momentum flux to the gas pressure. Then 
at the sonic radius S s = eyc/c s (r s ) (see straight lines in Fig . [Tb . 
In a chemical equilibrium determined purely by two-body col- 
lisional processes and photoionization, the equilibrium abun- 
dances Xi = n i/ n °f a U species are functions of the tempera- 
ture and F v /n only, where F v is the radiation flux at frequency 
v such that J F„dv = F =L/Anr 2 . On the other hand, the tem- 
perature r eq arising from the equilibrium of photoheating and 
two-body collisional cooling is a function of Xi an d F v /n only. 
Therefore, for a particular SED /„ = F v /F, the equilibrium 
temperature is determined only by L/4irr 2 n, i.e., T sq lies in 
one-to-one relation with 5 and we can write r eq = 7^(5;/,,). 
For particularly hard spectra, instead of being one-to-one, the 
function r eq (S;/„) can be multivalued in a certain range of S. 
Then, the temperature of a fluid element depends on its ther- 
mal history. Since thermodynamic perturbations at constant 5 
are isobaric, the possibly multivalued function r eq (S) at fixed 
/„ determines the thermal phase structure. 

For a fixed SED, r eq (S) is sensitive to the metallicity 
for 10 5 K < T < 10 6 K because helium line cooling and 
Bremsstrahlung cooling dominate the cooling rate at any 
metallicity (with a comparable contribution from iron at high 
temperatures in the metal-enriched case), while oxygen and 
iron photoionization heating are by far the most important 
heating pro cesses in the metal rich case (e.g., Kallman & Mc- 
|Cray|l982j ). Figure [I] shows that for the /„ oc v~ l J spectrum 
for 0.1 Ryd < hv < 10 3 Ryd, at near-solar metallicities, the 
Compton-heated hot coronal phase appears at S ~ 10 (as is 
well known, see, e.g., Krolik et al.||1981 1, whereas in a gas 
with subsolar metallicity Z < 0.1 Zq, the coronal phase does 
not appear until the ionization parameter reaches 5 ~ 10 3 . 
The optically thin r eq (S) becomes independent of metallicity 
atZ~O.OlZ . 

In Figurelll we also show the relations H = ec(^ffim p /kT) l l 2 
that must hold at the sonic radius for three values of the radia- 
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FIG. 1 . — Metallicity dependence of the local thermal and statistical equi- 
librium temperature of a photoionized gas under optically thin conditions as 
a function of the ionization parameter 3 (see text). The functions 7^(3) for 
four gas metallicities, expressed in units of the solar metallicity, were cal- 
culated with the photoionization code XSTAR (Kallman & Bautista 2001 1 
for an f v oc u -1 spectrum between 0.1 Ryd and 1,000 Ryd. The curves 
r K] (H) are single valued and thus all equilibria are stable for the particular 
choice of spectrum, but they need not all be stable for harder spectra. The 
plot shows that 7^(3) is nearly independent of metallicity for Z < 0.01 Zq, 
and that a hot coronal phase appears at 3 ~ 10 3 in a gas with Z < 0.1 Zq 
and at 3 ~ 10 in a more metal rich gas. The straight lines are the relations 
3 = eciyfintp/kT) 1 ' 2 that must hold at the sonic radius for values of the ra- 
diative efficiency, from left to right, of e = 0.025, e = 0. 1, and e = 0.4. The 
temperature at the sonic radius is found at the intersection with the r e q(S) 
curve. 



tive efficiency: e = 0.025, e = 0.1, and e = 0.4, which, spec- 
ulatively, might be expected for a rapidly rotating black hole 
with a geometrically thin retrograde disk, a nonrotating black 
hole with a thin disk, and a rapidly rotating black hole with 
a geometrically thin prograde disk, respectively (see, e.g., 
Novikov & Thorne 1973, Zhang et al.| 
[20061 120081 INoble et al.|2008| and refi 



I Novikov & Thorne|1973||Zhang et al.|1997||Beckwith et al. 
| 2006| |2008| |Noble et al.|2008| and references therein). The 
temperature at the sonic radius is found at the intersection 
with the r eq (S) curve, that is, we have the implicit relation 

T s = r eq [S s ;/„(r s )] = T e<i [ec{^^m v /kT s ) l l 2 ;f v ], which can be 
solved for the temperature at the sonic radius T s as a function 
of e and /„ (e.g., |Ostriker et al.||T976| >. Evidently, in metal 
poor gas, the equilibrium temperature at the sonic radius is a 
strong function of the radiative efficiency; highly radiatively 
efficient accretion is susceptible to thermal runaway, where 
the ionization state converges to full ionization as the Comp- 
ton heating overtakes t hermal evolu tion and the gas becomes 
fully ionized (see, e.g., Kroli k|1999| and references therein). 
Note, however, that the gas may not attain the Compton tem- 
perature if the inflow time becomes shorter than the heating 
time for combined Compton and photoionization heating (see 
§ |2.2| i. Given T s and the sonic radius determined from the rela- 
tion involving c s (r s ) and the forces acting on the gas, we have 
that r s = (l-£)fim p GM Bli /2'ykT s . For example, for power law 
SEDs F v oc v~ l - 5 and radiative efficiencies e ~ 0.1, typical 
temperatures at the sonic radius are T s ~ 10 5 K, and so the 
sonic radii are r s ~3x 10 14 M2 cm. 

We can write the usual ion ization parameter £ = L/r 2 n = 
AirkTcE of Tarter et al. ( 1969 1 2 as £ = 4irvnm p L/M, which in 
the optically thin limit becomes £ = 4Trv/j,m p ec 2 . In this limit, 

2 |Tarter et al.| jl969| define § in terms of the hydrogen density «h and 
we do so in terms of the total particle density n; in what follows, we do not 
always rigorously track the dependence on the molecular weight fj,. 



the temperature T eq , which is in one-to-one relation with £, is 
also in one-to-one relation with the velocity, r eq (v;/„,e), and 
so one can write equation Q as a differential equation with 
a single unknown function v(r). When the flow is isother- 
mal, we have the well-known asymptotic solution far from 
the black hole 



,3/2 



— I c s (r s ) (r ^ r s , isothermal). (3) 



If the temperature of the photoionized flow can decrease with 
radius and thus the infalling gas can acquire momentum be- 
fore it heats to ~ T s , the asymptotic velocity v(r) far from the 
black hole can exceed the isothermal value given in equation 
([3]l by a factor of several. We thus write 

4.5 T ( - ! c s (r s ) (r»r s ), (4) 



where T > 1 . 

Ideally, we would like to match this solution to the condi- 
tions far from the black hole, where the density and the total 
pressure are n^ and P^. One could attempt to set the bound- 
ary conditions n(r x ) = n^ and n(r x )kT eq {rx) = P^ = n^kT^ 
at some radius r\, where the last relation defines T^. This 
may indeed be possible at very low densities or very low ra- 
diative efficiencies. At high densities and efficiencies, how- 
ever, the photoionization equilibrium at some well defined 
radius rj on ^> r s abruptly transitions into a neutral state, i.e., 
on its way toward the black hole the gas passes a quasista- 
tionary ionization front. Gas density, velocity, temperature, 
and pressure may be discontinuous at the ionization front; 
furthermore, the warm or cold neutral gas in the immediate 
vicinity of the ionized region may be supersonically turbu- 
lent, in which case the gas density is inhomogeneous and the 
ram pressure of turbulent flows cannot be neglected. We pro- 
ceed to an attempt to determine under which conditions is a 
steady, strictly time-independent accretion across the station- 
ary ionization front possible. 

2.2. Size of the H II Region around the Black Hole 

We will assume that the ionized region is surrounded by 
warm, partially-ionized gas, though in reality, the accreting 
black hole and its ionization sphere may be embedded in a 
cold, molecular, and supersonically turbulent medium. How- 
ever, since the UV radiation from the black hole and the 
He II recombination radiation from the photoionization annu- 
lus (see § |4.2| below) will dissociate molecules in the vicinity 
of the H II region, a layer of warm, atomic gas should sur- 
round the ionization sphere even if the molecular phase exists 
at somewhat larger optical depths. We further assume that far 
from the sonic radius, where the black hole's gravity can be 
ignored, the ionized gas is in gas pressure equilibrium with 
the surroundings, i.e., 



T m "mi = /turb 7k «m (f > r s ), 



(5) 



where «hii is the total density of ions and electrons in the H II 
region, n^i is the atomic density in the neutral gas, Tun and 
Tm are the respective temperatures, and / tur b > 1 is a factor 
quantifying the degree of pressure enhancement due to turbu- 
lence in the neutral gas. Pressure equilibrium will be violated 
in non-steady-state, episodic accretion (§|3]l. 

Depending on the SED of the central source, the photoion- 
ization rate will be dominated by primary photoionizations or 
by secondary photoionizations carried out by photoelectrons 
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(e.g.,|Shull & van Steenberg|1985[|Xu & McCray|199lT|DaF 
|garno et al.||1999| l. Taking into account only the ionization 
of hydrogen from the ground state, the total rate of photoion- 
ization in the annulus can be written N[ on = fi on L/Eu, where 
E H = 13.6eV,L = 10 40 L 40 erg s 1 is the luminosity shortward 
of En, and /;„„ ~ (Eh/ hv) is the average fraction of the energy 
of an absorbed photon that goes into photoionization. For an 
almost fully ionized gas and a power law spectrum f u oc v~ a , 
this fraction equals /i on = (a— 1) /a for a > 1; in what follows, 
it should be borne in mind that /; on depends on the shape of 
the SED. 

In a photoionization equilibrium N[ on equals the total 
hydrogen recombination rate in the ionized gas N rec ~ 

f 7 r 'fon a «( 7 Hn)"H+«e, where a B (T) — 2.6 x 10" 13 r 4 _1 cm 3 s" 1 
is the approximate hydrogen recombination coefficient in 
the on -the-spot approxima tion at te mperatures 10 4 K < T < 
10 5 K ( |Ferlandet al.|l992| quoted in |Jappsen et al.|2007| , and 
«h+ ~ n e ~ jnim are the ion and electron densities, respec- 
tively. Though not entirely justified in Strtimgren spheres, 
the on-the-spot approximation is reasonable unless the den- 
sity profile within the ionized region is sharply peaked toward 
the center (Ritzerveld 2005}, as is assumed not to be the case 
here. 

Equating the ionization rate to the recombination rate we 
obtain 



4.4x 10 1 



,1/3 ,1/3 
J ion ^40 



HII.4.7 



,2/3 2/3 -2/3 



(6) 




HI.3.7 



atomic 

s 

^'molecules? 
turbulence? 



FIG. 2. — Schematic representation (not to scale) of the structure of the 
H II region surrounding an accreting seed black hole. For normal proto- 
galactic densities, the radius of the H II region r lon (eq. ^J) is a few orders of 
magnitude larger than the sonic radius r s (eq. fTI ). Ga s is m local thermal and 
statistical equilibrium outside the radius r c q u j (§|2.3|, which may be slightly 
larger or smalle r tha n r s . The ionized gas contains a small, nonzero neu- 
tral fraction (eq. |20| ), which gives rise to photoionization radiation pressure. 
The ionized region is surrounded by a shell of at omic gas of thickness Ar,ji ss 
where molecule photodissociation is efficient (§ |4.2| . The outer, molecular 
shell may be supersonically turbulent. 



where T Hll = 5 x 10 4 7^,4.7 K, T m = 5 x 10 3 T HL3 j K, and 
"hi = 10 5 «5 cm" 3 (we justify our choice of the reference den- 



sity in § 2.4 below). If we express the luminosity in terms of 
the dimensionless ratio £ = L/ (4irGMBum p c j aj) of the lumi- 
nosity to the Eddington luminosity (here and henceforth, for 
Thomson scattering), the radius of the H II region becomes 

r „ 4 7 x 10 18 £ ' /3 fi£ Tml ^ J cm m 
r lon ~4./xiu /3 2/3 2/3 cm. (/) 

Vturb W 5 7 HI,3.7 

The ionization radius is normally much larger than the sonic 
radius, 

f\/3 A/3 j 



' ion 



2 x 10 4 



,•2/3 2/3 —2/3 
/turb'V r > 



iT St5 



(8) 



In Figure|2]we provide a schematic illustration of the structure 
of the H II region surrounding an accreting seed black hole. 

Self-gravity of the gas inside the H II region is negligible 
compared to that of the black hole, ^TrPrin+mp <C Mbh, when 
I < 0.054 ^/t ur b«57Hi,3.77Hn,4.7( r / r ion) _3 which is satisfied 
for r <C r loa , but need not always be true at the very edge 
of the H II region. There, however, the combined gravity of 
the black hole and the gas contained within r; on are negligi- 
ble compared to gas pressure gradients; the ratio of the Jeans 
length inside the H II region to its radius is normally much 

larger than unity, Aj/r ion ~ 28 f^t^nfT^^ 3 , 
and thus self-gravity of the ionized gas can be ignored. 

2.3. Conditions at the Sonic Radius 

We will assume throughout that the gas is almost fully ion- 
ized at the sonic radius, i.e., that r s <C r{ oa . This condition 
places an upper limit on the density of the neutral gas just 



outside the H II region, n <C «„ 



, where 



= 2.7x 10 



f l/2 fl/2 T 3/2 r 3/2 
11 ./ion * J HII,4.7 J s,5 



ft 



turb JHI.3.7 



M 2 



(9) 



and we have assumed that 7hh. 4 .7 also represents the temper- 
ature of the ionized gas at the sonic radius. This shows that 
for an isotropic central radiation source, the gas at the sonic 
radius is guaranteed to be ionized unless the accreting gas has 
densities far in excess of those expected for the diffuse proto- 
galactic medium and instead characteristic of self-gravitating, 
star-forming cores. 

The central luminosity can be related to the mass accretion 
rate by substituting equation (Q) into i = e^v/incaj / GM^h to 
obtain the dimensionless luminosity for steady state accretion 
through the H II region 



'0.001 



£-1 /turb Y Ml "5 7ffl,3.7 



fkii.4.7 T s5 



3/2 



(10) 



where e = 0.1 e_i is the radiative efficiency. Evidently, the 
mass accretion can be strongly suppressed by the heat- 
ing of the ionized gas near the sonic radius. The accre- 
tion rate is much smaller than the Bondi accretion rate 
A^Bondi = e 3//2 7rG 2 Mg H m ; ,«Hi/c 3 HI that would be calculated ig- 
noring photoionization and photoheating altogether, i.e., for 
radiatively-inefficient accretion 



'5.2 



e_i M 2 n 5 



.3/2 r 3/2 
turb I H13J 



(no photoheating/ionization) , (11) 



where have assumed an isothermal equation of state. 
To justify our assumption in § 2.1 that the photoionized gas 



is in local thermal and statistical equilibrium, we calculate 
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the ratio of the Bremsstrahlung cooling time t% 
10 U T 
radius 



- 2.5 x 

10 I1 r 1 / 2 «H I sat the sonic radius to the inflow time at the sonic 



to obtain 



c s (r s ) 



77 



1.45 ^yr, 



3/2 



s.5 



's.5 ^HIL4.7 



/turb T r H I,3.7 «5 M 2 



(12) 



(13) 



Cooling due to the recombination of H II and He III is 
comparable to and only slightly stronger than that due to 
Bremsstrahlung at T s < 2 x 10 s K. 

We also estimate the ratio of the Compton heating time 
at the sonic radius tc ~ 0.0675 GM fx 2 m e m p c / ££ 2 r s 7c, where 
Tc ~ 10 7 Tcj K is the Compton temperature, to the inflow time 
(see also Fig. 2 in |Sazonov et al.|2005) , 

7-1/2 



tc 



0.01 



i Te- 



rn 



At the relatively high gas densities considered here, which are 
required for a rapid growth of seed massive black holes in 
protogalaxies, the Compton cooling of the photoionized gas 



by the cosmic microwave background photons (see, 
cotti et al.|2008| > can be ignored. 



e.g.- 



Ri- 



Substituting equation ( 10 1 in equation ( fl"4| ) we obtain for 



ratio of the Compton heating time to the infall time 



11 



^ S 5 7ffll,4.7 



(15) 



£-1 /turb T «5 M 2 Tqj 7ffll,3.7 

13 1 and ( fl"5j ) suggest that heating and cooling 



Equations 

times at the sonic radius can be longer than the inflow time 
and that the ionized gas may not be in local thermal and statis- 
tical equilibrium at all radii r ^S> r s , especially if the accretion 
occurs below the Eddington-limited rate. However, the ion- 
ized gas should be in equilibrium at only a somewhat larger 
radius because the infall time increases rapidly with radius, 
rjv oc r 3 . 

So far we have ignored the dependence of the temperature 
at the sonic radius on the radiative efficiency, which as Figure 
[JJshows, can be strong. The precise form of the function T s (e), 
defined by the intersection of the T eq (S) curve and the T S (S) 
line, depends sensitively on the metallicity of the gas and on 
the SED of the central source and we do not attempt to model 
it in general. If the dependence can be approximated with a 
power low in a range of efficiencies, T s oc e , from equations 
( fD) , ( [Tol l, and we obtam 



^Brems 28 
OC 6 , 



oc e 



1-35/2 



(16) 



If, e.g., S ~ 1 for e > 0.1, we find that with an increasing ef- 
ficiency it becomes more difficult for the gas to achieve local 
thermal equilibrium at the sonic radius, but if the equilibrium 
is achieved, the central luminosity, and especially the accre- 
tion rate that is proportional to M oc l/e oc e -35 / 2 , decrease 
with increasing efficiency 

Ciot ti & OstrTker] 0Wf\ [200T1 12007) and |Sazonov et aT] 
(2005) have studied spherically-symmetric accretion of hot 
interstellar medium onto an X-ray quasar in an elliptical 
galaxy and have identified a limit cycle driven by Compton 
heating. The quasar heats the interstellar medium to tempera- 
tures exceeding the virial temperature, which leads to an out- 
flow and quenching of central accretion. This model differs 



from ours in that a metal-enriched environment is assumed, 
so that photoionization equilibrium temperature reaches the 
Compton temperature already at 5 ~ 50, whereas in our 
metal-poor model the transition to the Compton tempera- 
ture occurs at higher values of the photoionization para meter, 
5 ~ 10 3 . Therefore, in the model of |Sazonov et al.| ( |2005| ), the 
ionized gas can reach the Compton temperature well outside 
the sonic radius. Another crucial difference is the assumed 
ionized gas density profile far from the sonic radius: the hot 
gas surrounding a quasar was assumed to be hydrostatically 
confined by the gravity of the host galactic stellar spheroid 
such that its density declines steeply with radius, «hh oc r~ 2 , 
which implies that H is roughly independent of radius and so 
the Compton-heated equilibrium can exist at arbitrarily large 
radii. In our model, since the H II region surrounding a seed 
black hole is confined by external pressure at radii /■»/■; and 
is roughly isothermal (consistent with the very weak depen- 
dence of r eq on the dimensionless ionization parameter for 
10 ^ S < 10 3 ; see Fig. [I]), the density inside it is approx- 
imately independent of radius. Then 5 oc r" 2 , and the gas 
is progressively farther from being able to heat the Compton 
temperature at radii much larger than the sonic radius. 



2.4. Protogalactic Density at the Edge of the H II Region 

We will now pause our investigation of radiative feedback 
effects to clarify our choice of the density of the ambient neu- 
tral gas surrounding the H II region. Rapidly growing proto- 
galaxies contain gas with a wide range of densities, and thus 
care must be taken to appropriately specify the gas density on 
scales relevant for regulation of accretion onto a seed black 
hole by radiative feedback effects. In the simplest picture, 
central gravitational collapse of the gas at the center of a pro- 
togalaxy produces a distribution in which gas density is spher- 
ically symmetric and a function of radius only, n(r). If the gas 
is approximately isothermal and quasi-hydrostatic in the pres- 
ence of turbulence, n(r) oc r~ 2 . If we ignore the gravity of any 
seed black hole and let / turr) denote a radius-independent tur- 
bulent pressure enhancement, the density scales as 



n(r)> 



/turb C 2 

27r G m p r 2 
'6x 10 4 /t^ Tin, 



3.7 



r 

1 pc 



cm 



(17) 



Equation 
profile resulting 



17 1 is a good approximation to the density "cusp" 
from central gravitational collapse in a 
10 8 M m cosmological halo in the simulations of 



central 

Mhaio ~ IO^Mq cosmological halo in 
Bromm & Loeb (2003) and in the simulations of Wis e et aL| 
( |2008| ), who found turbulent Mach numbers M ~ 3 in the 
cusp, which would imply / tur b ~ 10. 

If a seed black hole is located near the center of the density 
cusp, radiative effects may prevent the gravitational collapse 
from proceeding to the arbitrarily large densities. We will here 
assume that radiative effects prevent collapse at radii smaller 
than the radius of the H II region, i.e., that the density pro- 
file of the neutral gas has a "core" on scales ~ r lon such that 
density within the H II is roughly uniform far from the sonic 
radius. For self-consistency, we substitute n = n(rj on ) from 
equation ( 17 1 in equation (j7]) and solve for rj on , and in turn for 



8 



ACCRETION ONTO SEED BLACK HOLES 



«(>"ionX to obtain 



i tu 



o q i rvl7 •'turb HI, 3. 7 
' lon.cusp J J.O X 1U — YtT> a/T 
/ion *HII4.7 2 



-4.5 xlO 6 ' /,on 



2 T 6 M 2 
^11.4.7 ™2 



f T 1 



(18) 



'turb J HI.3.7 

These estimates are but crude self-consistency conditions and 
suffer from a strong sensitivity to the turbulent Mach number 
and other parameters. The estimate of gas density at the edge 
of the H II region given in equation ( fT~8] > can only be used as a 
rough guide for the range of densities that should be addressed 
in the ensuing analysis. Our choice of the reference density, 
10 5 cm" 3 , is compatible with the self-consistency conditions 
in a protogalaxy in a cosmological halo of Mhaio ~ 10 8 M Q for 
/ion ~ 5 , ignoring any recent supernova activity in the center 
of the halo, which can drastically reduce the central density 
(see, e.g., Wada & V enkatesan 2003 Kitayama & Yoshida 
2005] | Wise & Abel|2008| >. " 

Equations (|18[> also suggest that "minihalos" with masses 
Mhaio ~ 10 6 M Q and Mach numbers M > 1 will be fully ion- 
ized out to r > 100 pc if, somehow, r? max ~ 1 is realized at an 
early instant prior to the expansion of the H II region, and the 
interior of the H II region has not had chance to heat beyond 
10 4 K. This is seen in the simulations of radiative feedback 
during accretion onto seed black holes in cosmological mini- 
halos of |Alvarez et aE1p008l >. 

2.5. Continuum Radiation Pressure 

Here we estimate the effects of the continuum radiation 
pressure in the interior of the H II region. The continuum radi- 
ation pressure acceleration is a ra( j = <7H£xHi/47rr 2 m p c, where 
Xhi is the abundance of neutral hydrogen in the ionized gas, 
and (7h = F~ l j a^F^dv is the frequency-averaged mean ab- 
sorption cross section. In a highly ionized gas, most of the 
photoelectron energy goes into heating, so in ionization bal- 
ance 

Qfi(?Hii) 

L o- h /(4tt r 2 n e E H ) + a ion (T m ) ' 

where <r H = F~ l Eu f (hu)~ l a^ u F u dv and ai on (T) is the colli- 
sional ionization rate. At temperatures <C 10 5 K, photoion- 
ization typically dominates collisional ionization. The neutral 
hydrogen abundance then becomes 

,2/3 / x 2 



XHI ~ 7.6 X 10" 



ft 



on f. 



1/3 M 1/3 « 1/3 

turb m 2 n 5 



r l/3 
HI, 3. 7 



r 

Hon 



(20) 



which implies that at radii not much smaller than r mn , the 
neutral abundance is > 10 -6 , and in this regime, since cjh ~ 
10 6 e>T, radiation pressure due to photoionization exceeds that 
due to Thomson scattering. The acceleration due to the for- 
mer is then given by 

oh a B (Tm) n e E H 

a m d=~ • (21) 

<j H m p c 

Note that, unlike in the case of Thomson scattering, this is 
independent of distance from the black hole because here the 
abundance of absorbers increases with the square of the ra- 
dius, xhi oc r 2 . The cross section ratio is cth/oh ~ (3 + a)/(2+ 
a) for power-law spectra F v oc v~ a . We compare a rac j to the 
black hole's gravity a grav = -GM^/r 2 to find, 



where we have assumed a fiducial a = 1 .5 spectrum, implying 
/ion ~ 3 and (Th/cth ~ %■ The relative importance of photoion- 
ization radiation pressure depends on distance, and dominates 
at larger radii. This shows that for a wide range of param- 
eters, for near-Eddington accretion, the radiation pressure in 
the H II region at radii r < r lon greatly exceeds the black hole's 
gravity, which implies that steady state accretion may occur at 
a reduced rate. Inward accretion is still possible if a positive 
pressure gradient is set up to counteract the radiation pressure 
force. 

To estimate the reduction of the steady-state accretion rate 
as a result of the photoionization radiation pressure in the 
H II region, we simplify the mathematical problem by ignor- 
ing any heating of the gas near the sonic radius and setting 
T s = 7hh = const; this simplification, as we shall see, will turn 
out to be unreasonably optimistic. Then, after a change of 
variables, 



w \ 



c s ,hii 



y : 



GMbh 

2 C s,HII 



r 



(23) 



equation ( 2} describing momentum conservation in the accre- 
tion flow becomes 



dxv 
dy 



(w — 1) = w 



2 2 



(24) 



j y 

where <p is a dimensionless parameter proportional to the pho- 
toionization radiation pressure acceleration, 

, GM BH 
-T— "rad 



2c! 



'2x 10 



-5 /turb «5 M 2 7hii,3.-, 



(25) 



Equation ( |24| >, subject to the regularity condition imposed at 
the sonic radius, w(y s ) = 1, where y s = 4>~ x [(\ +2<p) l l 2 - 1] ps 1, 
can be solved in closed form in terms of the transcendental 
Lambert W(x) function, defined implicitly via x = We w , to ob- 
tain 



w(y) = 



\ 



-W 



exp 



2(# 2 + 2) | 2(# 2 + 2) i 



iy exp | 



# 2 -3y/2 + 2 



y 



(26) 



where in the second line we applied the expansion W(x) = 
x+0(x 2 ). 

For <fi^il, the ratio of density n oc (r 2 v) _1 oc (j^vv) -1 at ra- 
dius r ^> r s to that at the sonic radius has asymptotic form 



n(r s ) 



r,<t> r/r s 



n(r) 
n(r s ) 



(27) 



0=0 



This result is exact for isothermal accretion, but here, we 
adopt it as an optimistic estimate of the density reduction due 
to photoionization radiation pressure in the case of realistic, 
nonisothermal accretion. 

We expect significant reduction in the accretion rate when 
"hii^s) "C HHii^ion) or, according to equation (27 i, when the 
quantity 



f l/3 A/3 „ 1/3 1/3 



T . 1/3 ,.1/3 
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is comparable to or greater than unity. The dimensionless ac- 
cretion rate can be eliminated by substituting equation ( [T0| 
that relates the accretion rate to the conditions at the sonic 
radius, to obtain 



0.04 



f l/3 f l/3 T i/ 3 1/3 r 2/3 1/2 2/3 . . 
/ion /turb 1 e -l y HI,3.7 i s,5 W 5 ™: 



2/3 
2 



r 10/3 
'mi,4.7 



(29) 

Thus we find that accretion will proceed in one of the two 
distinct regimes. When 



«1, 



(30) 



radiation pressure does not affect the accretion rate even if the 
combined gravitational and radiative acceleration is outward 
in parts of the H II region. However, when 



>1, 



(31) 



the ionized gas assumes a density gradient at radii r^s> r s such 
that the density increases outward, toward the edge of the H II 
region. In this case, central accretion is suppressed at least by 
the factor 

£ (x exp (-/-^j , (32) 

and probably more, because the central drop in density will 
lead to a photoionization equilibrium at a higher tempera- 
ture, thereby increasing gas temperature near the sonic radius. 
Therefore, for (j)r lon /r s > 1, radiation pressure acting on the 
ionized gas imposes a more restrictive limit on the accretion 
rate than the mere photoheating near the sonic radius. The 
two limits are complementary: while the photoheating limit is 
more restrictive for small black hole masses, ambient densi- 
ties, and turbulent Mach numbers, the radiation pressure limit 
is more restrictive in the very opposite regime. A combina- 
tion of these two limits suggests that a steady state solution 
for near-Eddington accretion does not seem to exist. 3 

A fraction / res < 1 of the luminosity emitted by the central 
source passes unabsorbed through the H II region and reaches 
its edge, where most of it is absorbed in the partially ionized 
gas. Momentum deposition by the residual luminosity leads 
to a density drop by the factor of e\p(-ip) where 



/res ^ 



4tt r ion 7 n H n k T ml c 

/res 7hII,4.7 \ ± ^ion 



0.033 



/ion T s i 



(33) 



Since T s > T m i, radiation pressure in the partially-ionized 
edge of the H II region is always less important than in its 
interior. 

2.6. Accretion at Low Radiative Efficiencies 

In this work we focus entirely on the regime in which cen- 
tral accretion is radiatively efficient, e.g., e 0.01, as it 
should be the case when the radiation is emitted by a geo- 
metrically thin disk. In this regime, heating of the gas at the 
sonic radius to > 10 5 K, and perhaps even to the Compton 
temperature at high metallicities, seems inevitable. We hope 
to generalize our analysis to radiatively inefficient accretion, 



3 |Wang et al.|{2006| concluded that radiation pressure was negligible be- 
cause they erroneously assumed complete ionization (xm = 0). 



which should be accompanied by outflows and shadowing, in 
a subsequent investigation. Here, we will make only a brief 
mention of this alternate regime. 

If we ignore the mechanical and radiative effects of the out- 
flows and the shadowing, which is an implausible and incon- 
sistent assumption, a change in the thermodynamic behavior 
occurs at e < 10~ 3 , namely, at such low efficiencies the gas 
at the sonic radius and throughout the H II re gion can remain 
at temperatures well below 10 5 K (see, e.g., |Buff & McCray 
T974| [Hatchett e t al.|1976HKrolik et al.|1981[|Kallman & Mc- 
Cray|il982||Krolik & Kallman|1984||Lepp et al.|1985||Don- 



ahue & Shuhj|l99l| l. If, e.£7T s ~ r Hn ~ 10* K _ and t^T; 
we find that, still, under a wide range of conditions character- 
istic of protogalactic clouds, an extended (ri on ^S> r s ) H II re- 
gion exists around the seed black hole. In this case, photoion- 
ization radiation pressure suppresses central accretion when 
4> r ion/r s > 1, which can be expressed as a condition on the 
radiative efficiency; the suppression takes place when e > e cr ; t 
where 

0.005 

£crit ~ 7 f 2 T 2 n^M 2 ' 

Jion 7turb J HI,3.7 "5 m 2 

For e > e cr it, the ionized gas may not be able to heat above 
10 4 K, but the photoionization radiation pressure suppresses 



accretion well below the rate £ s s . derived in equation ( 10 1. For 
e *C e C rit, accretion is limited by conditions at the sonic radius 
(§ |2.3| l. The accretion rate may become close to the "Bondi" 
rate calculated ignoring radiative effects altogether (eq. fTT) ) 
if the neutral gas surrounding the H II region is weakly super- 
sonically turbulent, but the ionized gas inside it is not. 

We would like to reiterate that the analysis presented in 
this subsection is somewhat unrealistic, because radiatively- 
inefficient accretion flows are probably accompanied by out- 
flows that can radically alter the hydrodynamic, thermo- 
dynamic, and chemical structure of the region surrounding 
a seed black hole. Also, the radiation field produced by 
a radiatively-inefficient accretion flow is potentially highly 
anisotropic. 

2.7. Lya Radiation Pressure 

Another source of pressure is from the Lya photons that are 
produced throughout the H II region and can become trapped 
within the ionized region and the surrounding neutral shell 
(see, e.g., |McKee & Tan"l|2008| a nd references therein). Let 



Osterbrock 



1989) denote the fraction of the 



/Lya ~ 3 (e.g., 

total power emitted by tne central source above E H that is 
reprocessed into Lya photons. Consider a shell of neutral 
gas at radii r; on < r < r[ on + Ar s h e ii, where Ar sne ii denotes the 
thickness of the shell. For shells with Ar sne ii < Hon, the line- 
center optical depth of the shell To = tJo(!r)«HAr s h e ii, where 

a (T) = 5.9 x 10- I4 r 4 " 1/2 cm 2 , is large 

TO ~4xio-°^ <1/3 <f WA^Y (35) 

\ Hon / 



.2/3 T 7/6 
turb 7 HI,3.7 



Since the shell is very optically thick, a photon injected at the 
edge (r = r; on ) near line center will scatter across the ionized 
region many times before escaping the shell. The ionized gas 
within the H II region has some residual optical depth due to 
the presence of a small neutral fraction. Under the conditions 
considered here, the ionized gas is optically thick to the pho- 
tons in the core of the Lya line, and is marginally optically 
thin to the photons in the wings. We ignore this complication 
and treat the ionized gas as optically thin. 
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To estimate the number of times a photon injected at the 
edge of the neutral gas crosses the H II region, we assume that 
Ar s heii < r[ on , and carry out a simple Monte Carlo calculation 
of a photon's frequency diffusion before it escapes the shell. 
For this, we employ the accurate expressions for the trans- 
mission coefficient and the reflection frequency redistribution 
function at the shell edge that Hansen & Oh| ( 2006 1 obtained 
from Monte-Carlo resonant line scattering calculations. In the 
dust-free limit, the number of times a photon injected near the 
line core reflects against the walls of the ionized region is ac- 
curately approximated with 



Reflect « 0.609 (a r ) U My e 



0.659 -0.00607 [ln(2120.0 a T )f 



(36) 



where a(T) = vi^jlv^^ = 4.72 x 10 4 T 4 is the ratio of 
the natural to the Doppler line wid t h. This is c lose to the 

Reflect ~ 



1975 1. 



asymptotic formula due to Adams ( 1972 
15 (tq/10 5,5 ) 1 / 3 , in the optically thick limit (quoted from|Dijk- 
20081. AtT 



stra & Loeb 



10 4 K,forr = 



(10 9 , 10 10 , 10 ir Vthe 
numbers of reflections take values N re n ect ~ (250,610, 1400). 
Rayleigh scattering contributes to the opacity negligibly for 
To < 10 10 . Destruction of Lya photons by stimulated two pho- 
ton emission is negligible in the range of optical depths that 
we consider. Destruction by H" ions in the partially ionized 
shell surrounding the H II region is potentially important (see 
§ |4.2| i, but because of significant uncertainties, we ignore it 
and adopt Adams' formula in what follows. 

If the radiation pressure drives a coherent expansion (out- 
flow) in the neutral shell, this may promote photon escape 



and reduce P\v a /Pt 



gradient of magnitude Av, Bonilha et al. ( 1979 



If the outflow contains a velocity 

see also 



Bithell 1990 1 estimate that /V sca tter is reduced by a factor 
~ (1+0.04 lAv/vDopI 3 / 2 )" 1 , where 

v'Dop ~ c s ~ 10 km s is 
the thermal or turbulent Doppler velocity. The same reduc- 
tion factor should apply to A^ re fl e ct- Thus, even for a highly 
supersonic outflow Av/vDop ~ 8, the number of reflections to 
escape and the Lya radiation pressure to which it is propor- 
tional will drop by only a half. 

The energy density in Lya photons can be approximated 
with (see, e.g., |Bithell|199^|Haelmelt|1995] l 



t/i 



Lya 



reflect 



47r(r+Ar she ii) 2 c 



(37) 



From this, the radiation pressure P]_y a = \U]_ ya scale height is 
~ j Hon- At depth equal to one scale height we have, with the 
fiducial choice f Lya = |, 



*Lya 
P 

Jgas 



'3.8 



4/9,4/9 



f 5/9 T l/18 r 5/3 
J io" 'win ' 1 1 



(38) 



HI, 3. 7 ^HII.4.7 



For near-Eddington accretion the Lya pressure can signifi- 
cantly exceed the gas pressure and impart an outward impulse 
to the neutral gas surrounding the H II region. This result ap- 
pears to imply that even at relatively small I, the pressure of 
the trapped resonance line radiation remains in excess of the 
gas pressure. 



If we speculatively require Plto 
accretion, we invert equation (pSl 
the accretion rate, I < 4rit,Lya, where 



< ^gas for strictly stationary 
to derive an upper limit on 



^crit.Lya 



:0.05 



,5/4 T l/8 



,15/4 



HI.3.7 ' t HII,4.7 



ftJrt M 2 "5 



(39) 



Because of the strong dependence on Tun and our not hav- 
ing sought, within the confines of the present work, a self- 
consistent dynamical solution for gas flow in the presence of 
resonance line scattering radiation, we are not able determine 
whether the continuum photoionization pressure or the res- 
onance line scattering is more constraining to the maximum 
accretion rate for strictly stationary quasiradial accretion that 
can be achieved. Additional sources of uncertainty affecting 
any attempt to estimate the impact of resonance line radia- 
tion pressure are the topology of the density field beyond the 
spherically symmetric approximation, as the radiation may 
"leak out" eve n through a relatively small hole in the neutral 
shell (see, e.g., |Hansen & Oh|2006| [McKee & Tan|2008| l and 
the dynamical stability of shells accelerated by resonance line 
radiation (Ma thews|1992| l. 

3. EPISODIC ACCRETION 

3.1. Mechanics of Inflow and Outflow 

Here we revoke the assumption of strictly stationary flow 
and attempt to describe a sequence of periodically recurring 
stages in the seed black hole's accretion cycle. Let £ max < 1 
denote the peak luminosity in units of the Eddington luminos- 
ity. The peak luminosity is reached due a sudden infall of ma- 
terial; the thermodynamic state of the infalling gas is assumed 
to be set by a central luminosity well below the peak luminos- 
ity. The infalling gas reacts to the sudden rise in central lumi- 
nosity on a finite time scale. Facing a sudden rise in the cen- 
tral luminosity, the accreting gas heats by photoionization and 
gets accelerated outward by the radiation pressure. The pho- 
toionization heating time f p hoto = 47rr 2 (7~ l)' 1 kTun/ ^xh&hL, 
where cth is the average photoionization cross section and \n 
is the neutral hydrogen abundance (see §|2.5[l, is given by 



fphoto °° 2.9 



r 3 

i HII4.7 



Jturb /epi,l «5 7ffl.3.7 



yr, 



(40) 



where we have included a potentially large density enhance- 
ment /-pi = 10 / e pi.i > 1, due to, e.g., infall or a lack of pressure 
equilibrium, over the steady-state density of the H II region. 
The photoionization heating time scale is thus rather short and 
comparable to the inflow time at r ~ r s . 

Following the sudden increase in central luminosity, mate- 
rial in the outer part of the H II region, r ~ r lon , is accelerated 
outward by the radiation pressure if a la d > ^ > gas.Hii/'"M m /> n Hii, 
which implies the condition (cf. eq. (28j) 



1/3 ,1/3 



3.3 



/ion ft 



(,1/3 ^1/3 „l/3 M l/3 



uirb /epi.l ^m'ax ^111,3.7 U 5 



T 2 

J Hn,4.7 



> 1. 



(41) 



Photoionization heating and radiation pressure will act to 
erase the infall-induced central density enhancement and 
drive the central density and the accretion rate down, toward 
the limits imposed by thermodynamics (§ |2.3| l and radiation 
pressure (§ |2.5| >. The accretion flow remains directed inward 
close to the sonic radius, but far from the sonic radius, the 
inflow gives way to an outflow. 

If the outflowing material reaches radii ~ r; on , it encoun- 
ters a shell of denser hot gas that has been photoevaporated 
from edge of the H II region. The hot shell is overpressured 
by a factor ~ 7m//turb7ffl and is expanding in both radial di- 
rections. As the central density and luminosity drop, radia- 
tion pressure is not able anymore to accelerate the ionized gas 
against the positive pressure gradient in the exterior of the H II 
region. Then the pressure gradient resulting from the rapid 
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decline of radiation pressure and from the photoevaporative 
heating accelerates the gas near the edge of the H II region in- 
ward, in an implosion back toward the black hole. The gravi- 
tational acceleration is subdominant compared to the pressure 
gradient acceleration, |a pres | > a ra d.max > |a gr av| (see eq. p2) ) 
until the infalling gas approaches to within a few sonic radii 
from the black hole, where, by definition, gravity becomes 
competitive with pressure. The returning gas acquires a ve- 
locity similar to the local sound speed of the photoevaporated 
gas, which is perhaps somewhat larger than e Si Hn('"ioa)> an d 
returns to the vicinity of the black hole on a time scale 



C S ,HII 



4x 10 4 



(,1/3 



HII.4.7 



,2/3 2/3 T 2/3 
/tuib w 5 7 HI,3.7 



yr. 



(42) 



The return completes the accretion cycle. The return time may 
be shorter by the factor of a few than this estimate if the over- 
pressuring during to photoevaporation, which we ignored in 
equation ( |4"2| ), is taken into account. 

The pressure of the Lya scattering radiation trapped within, 
and in the neighborhood, of the H II region may further reduce 
the average accretion rate by imparting an outward impulse to 
partially ionized gas and thus possibly significantly extending 
the return time f re tum- In view of the lingering concerns re- 
garding the transfer of Lya radiation that we have raised in 
§ |2.7| we ignore it in the following estimate of the duty cycle. 

3.2. Estimates of the Duty Cycle 

If, ultimately, after the central density and luminosity have 
decreased substantially, thermodynamics imposes the most 
stringent limit on the accretion rate, we might expect that the 
luminosity drop toward the steady-state rate ~ 4.s. derived in 
equation (TO I. Because of the sensitivity to the various pa- 
rameters of our model, we are not able to generally determine 
whether the luminosity reaches 4. s . in time f re tum separating 
consecutive infall episodes. For the simplicity of the remain- 
ing analysis of episodic accretion, we assume that the mini- 
mum luminosity is indeed ~ 4.s. an d that the decline is expo- 
nential, 

/ Q \ f Aretum 

i{t) = i m .Ji^\ . (43) 



Then, the time-average luminosity is given by 



ln(£ max /4.s.) 

and the duty cycle, which we define as /duty = (£ 2 ) /(^) 2 (see, 



e.g., |Ciotti & Ostriker|2001| l, is in the limit £ s s <§C ^max given 

by 



/duty 



(45) 



ln(^max/4.s.) 4nax 

The weak logarithmic dependence on the ratio of the maxi 
mum to the minimum luminosity is deceptive; while we gen 



erally expect a duty cycle in the range / 



duty 



0.2-1, the av- 



erage accretion rate depends on the peak rate. We can attempt 
to generalize equation ( 10 1 to model the peak rate by replac- 
ing r s with 7hh (to mimic an initial absence of heating at the 
sonic radius) and by including the density enhancement factor 
/ ep i to obtain 

„ n t e -l /turb /epi.l T M 2 «5 ?H1,3.7 

^HII,4.7 



which would imply a duty cycle of /j uty ~ 0.6 (for / ep ; ~ 10) 
and an average accretion rate that is a factor of ~ 9 times 
higher than the steady state rate. 

The duty cycle estima t ed he re differs from the one 
derived by 



Ricotti et al. 



= (re/non) 1 / 3 , 

where for a stationary black hole r B = GM Bli /c 2 Hl ~2.4x 



(20081, / 



/duty .ROM 



10 16 Mi 7^/3 7 cm is the Bondi radius ignoring radiative feed- 



back ( Bondi & Hoyle|1944| l. This estimate assumes that the 
accretion is efficient only when r{ on <r B . In our picture, r B 
does not have a unique physical meaning because the accre- 
tion flow structure is strongly modified by photoionization, 
radiative heating, and the radiation pressure; episodic accre- 
tion through the H II region can proceed even when r; on ^> r B ■ 
The peak accretion rate £ max and the average rate (£) during 
episodic accretion can exceed the limits imposed by photo- 
heating and radiation pressure because of the presence of the 
finite inertia of dense infalling shells and the lack of pressure 
equilibrium in the gas photoionized from the inner edge of 
the H II region during luminosity maxima. We are not able to 
determine with certainty whether the accretion from a weakly 
turbulent quasiuniform density medium will be episodic, in 



which case it will proceed at the rate given by equations (44 



tion 



and (46 1, or whether it will be steady at the rate given in equa- 
T0| ) subject to radiation pressure-suppression derived in 



We proceed to discuss the implications of inhomogene- 



§[23 

ity and turbulence in the environment of a seed black hole. 



4. MOLECULAR COOLING AND CLUMPY ACCRETION 

4.1. Accretion of Self-Shielding Clumps 

As a protogalaxy grows, the baryonic inflow velocity into 
its center increases. Cosmological hydrodynamic simulations 
show that baryons accreting along the filaments of the cos- 
mic web can remain cold until they reach the central region 
of the protogalaxy. Supersonic baryonic inflow, in proto- 
galaxies with virial velocities > 10 km s -1 , drives turbulence 
in the galaxy. The resulting central turbulent Mach num- 
bers measured in the simulations of ~ 10 8 M cosmologi- 
cal halos are - 3 ( fWise et al.||2008l |Greif et aJ.||2008> . Su- 
personic turbulence is expected in view of the short cooling 
time f s hock.cooi ~ 10 rcj 1 yr of the gas heated at the termina- 
tion shocks of the cold infl ows (for somewhat higher veloc- 

Shapiro & Kang||l987| |Kang & 



ity, 50 km s shocks, e .g 

|Shapiro|[T992l see also |Gnat & Sternberg||2q08[ for depen 
dence on metallicity, which is weak for Z < 0.1 Z ). If 
the shocked gas cools on a dynamical time to temperatures 
~ 10 4 K, molecules form and the gas quickly cools further. 
Therefore, a fraction of the dense gas mass that is collecting 
at the center of the protogalaxy may reside at temperatures 
substantially below the 7hi ~ 5,000 Kthat we have somewhat 
arbitrarily taken as the temperature floor thus far. 

Near an accreting seed black hole, molecule formation 
could be enhanced even in the absence of dust. Hard X- 
rays emitted by the black hole could maintain a high elec- 
tron fraction in the neutral gas surrounding the H II region, 
which would catalyze molecule formation. Molecules can 
also form if traces of metals and dust are present. The cooling 
to temperatures characteristic of the molecular phase leads to 
an increase of turbulent Mach numbers. In supersonic tur- 
bulence, the local density exhibits intermittent strong fluctua- 
tions around the mean. Simulations of isothermal turbulence 
have shown that the density probability density function is a 
normal distribution in In p with mean (In p) = - 1 a 2 and disper- 
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sion a 2 = ln(l +b 2 M 2 ), where M. is the turbulent Machnum- 
ber and b w 0.26 for unmagnetized turbulence (e.g., Kritsuk 



et al. 2007| and references therein). For example, for M. = 10 
1% of the volume contains densities in excess of the aver- 
age density by a factor of 10 or greater. Dense clumps can 
enter the ionized sphere and remain self-shielded from pho- 
toionization and photodissociation. Since from equation d6]i 
the critical radius for photoionization is proportional to n" 2 ^\ 
a clump with an overdensity of 10 can withstand photoioniza- 
tion to radii ~ 0.2 r; on . 

The presence of neutral clumps in the ionized sphere may 
fundamentally alter the structure of the accretion flow if they 
are dense enough to survive photoionization and photoevap- 
oration, which requires densities violating the condition in 
equation ([9]). Such dense clumps may form as a product of 
turbulent fragmentation and may be self-gravitating and on 
their way to turn into star-forming cores, provided that they 
are not tidally disrupted by and accreted onto the black hole 
(e.g.,[Bonnell & Ric e|2008| . During accretion minima when 
I <C £ max , dense clumps may withstand photoionization and 
photoevaporation even at r < r s . In this case, they could 
accrete directly into the small optically thick accretion disk 
around the black hole. If the clumps occupy a substantial solid 
angle as seen from the black hole, the gas in their shadows re- 
combines, cools, and accelerates toward the black hole. 

We thus speculate that with the increase of turbulent Mach 
numbers and the progression of ever higher degree of clump- 
ing in the gas, the severity of the radiative feedback discussed 
in §[2]and §[3]decreases. The feedback-limited accretion may 
proceed at only a small fraction of the rate corresponding to 
the Eddington limit until turbulent inhomogeneities reach the 
level at which the densest clumps are self-shielded at ~ r s 
even at £ max ~ 1. Then, unless the gas supply is depleted by 
star formation and supernovae, the seed black hole may be 
able to grow efficiently and double its mass on the Salpeter 
time scale fsaip ~5x 10 7 e_i yr. Central stellar velocity dis- 
persions found in massive black hole-hostin g stellar systems 
in the local universe are > 30 km s -1 (e.g., |Barth et al.|2005[ 
and references therein); if massive black holes are not found 
in systems with smaller dispersions, this may be interpreted as 
a hint that efficient accretion commences only when the bary- 
onic gravitational potential well depth around the black hole 
exceeds a critical minimum value. 

4.2. Photodissociation Sphere around the Black Hole 

The existence of dense clumps in the neighborhood of the 
H II region may be contingent on the presence of molecules. 
The radiation from the black hole likely contains a photo- 
dissociating component below the Lyman edge. However, for 
a central source with a hard spectrum, two photon emission 
from the 2 1 S — > 1 1 S transition in the recombination of He II 
is a guaranteed source of H2 -dissociating photons. Photons 
in the Lyman- Werner (LW) bands of hydrogen are produced 
at the rate /Vlw 



5./ion,He£/£He (Johnson & Bromm 2007), 



where /ion , He is the fraction of p hoton energy that go ' elmfe 
helium ionization in mostly neutral gas, and £He = 24.6 eV. 
Some LW radiation may also be produced by the central 
source itself (Kuhlen & Madau 2005). If the gas surround- 
ing the black hole is dust-free and molecule synthesis is cat- 
alyzed primarily by H", one can estimate the dissociation 
depth Ardi ss in the neutral gas via 



47T 



where r<ji ss = r lon + Ardi ss , k^ is the rate for the reaction 
H~ + H — > H2+e~, and «h- is the equilibrium abundance of 
the H" ions. Adopting the temperature T = 1,000 K we h ave 



ku 2 « 1.2 x 10 cm 3 s" 1 and n H - ~7x 10" n e (see, e.g., Oh 
|& Haiman|2002| and references therein). In writing equation 
(47 1, we have ignored the potential enhancement of H2 for- 



mation rate in a supersonically turbulent gas, where molecule 
formation is particularly efficient in intermittent overdensities 
(see, e.g., |Pavlovski et al.|2002||Glover & Mac Low|2007] l. 

The maximum ionization fraction \e = n e /«H allowed if H2 
is to be dissociated in a sphere of radius 2r lon is obtained by 
substituting Ardi ss = r; on in equation ( |47] > and solving for \ e , 
to obtain the condition 



Xc 



<2.4x 10 



f 2 

-4 Jturb 



T 2 

i HI,3.7 



/ion T\ 



(48) 



HII.4.7 



where we have set /i on> He = 0.06. If the SED of the black 
hole extends into the hard X-rays, an electron fraction violat- 
ing the condition in equation ( |48| may be maintained in the 
neutral shell surrounding the H II region. This suggests that, 
even in an environment with primordial composition, the H II 
region associated with a seed black hole may be surrounded 
by only a thin (Ard; ss -C r; on ) photodissociation region where 
gas cooling and clumping is reduced (see also Johnson et al. 
2007 \ . Photodissociation may be even less important if the 



protogalaxy has been enriched with trace quantities of metals 
and dust (e.g., Omukai et al. 2008 ), as gas cooling can pro- 
ceed even in the absence of molecular hydrogen in that case. 
Also, if the accretion is episodic (§[3}, molecular gas can form 
during the quiescent periods in which central accretion is tem- 
porarily diminished or suspended by the radiative feedback. 

We have ignored the photodissociation of H", which re- 
duces the H2 formation rate by a factor S~ l , where S = 
1 +^7,H-CH-/47rr 2 A:H 2 nH, Nj,h- is the total dissociating pho- 
ton number luminosity, and oh- is the average photodisso- 
ciation cross section (see, e.g., Glover et al. 2006 Glover 
[2007l|Yoshida et al.|2007||Chuzhoy et al.|2007| and references 
therein). In our estimate of 7\Lh_, we will ignore the multiply- 
ing effect of Lya trapping (§ 2.7 1 on the H"-photodissociating 
photon number density; a more accurate approach, which 
would be too involved to include here, would be to solve 
for Lya resonance line and X-ray continuum transfer in the 
H" photodissociation region self-consistently. Chuzhoy et al. 
( |2007| ) estimate that for reprocessed ionizing radiation, the av- 
erage cross section per recombination photon times the aver- 
age number of photons per recombination varies in the range 
(an-) ~ (1.6-3.4) x 10" 17 cm 2 . Setting 7V 7 , H - ~ N KC , where 

^Vrec ~ (47r/3)aB(T[ni)«H + «e is the hydrogen recombination 
rate inside the H II region, we find 



5- 1 + 0.43 



,1/3 A/3 -4/3 M l/3 1/3 4/3 
1 /ion /turb M 2 n 5 i HI,3.7 



1 HII.4.7 



, (49) 



(47) 



which shows that H2 formation suppression is marginally 
important in the case of an accretion from a high density 
cloud. The suppression widens the photodissociation region 
surrounding the H II region. 

Having raised the possibility of a multiplying effect of Lya 
trapping on the photodissociation of H", we remark that the 
destruction of Lya photons by H", which is a processes that 
have ignored in § |2.7| may have to be considered in a self- 
consistent study of Lya transfer in a gas with primordial com- 
position. Because of significant uncertainties in the ionization 



MILOSAVLJEVIC, ET AL. 



13 



fraction and the H~ abundance in the neutral shell, which both 
depend on the SED are not calculated accurately in our sim- 
plified model, we attempt only a crude estimate. According 
to equation (7) of |Hansen & O h ( 2006 ), Lya photons scatter- 
ing through a shell of column depth N are destroyed before 
escaping if 



N>6Jx 10 



19 j hi 



,1/4 

HI.3.7 / &Lya,H- \ 

~JI* UO" 21 cm 2 / 

AH- 



-3/4 



(50) 



where CLya.H- ~ 4.4 x 10~ 18 cm" 2 is the H~ photodissociation 
cross section by Lya, and xh- is the H" abundance. Setting 
/V ~ n H Ar we find that Lya photons are destroyed before es- 
caping the shell of width Ar when the ionization fraction in 
the shell is 



X e > 0.0024 



„ ,8/9 T n/9 

d /turb 1 HI.3.7 



M/9 f*/9 4/9 M 4/9 ^4/3 
* /ion n 5 M 2 i HII,4.7 



Ar 



-4/3 



■(51) 



This crude estimate suggests that for accretion in a dense, 
dust-free medium with a central X-ray source that maintains 
a high ionization fraction in the neutral gas, Lya photons may 
be destroyed by H~ absorption before penetrating the shell. 
The issue warrants a self-consistent treatment of resonance 
line scattering and continuum radiation transfer coupled to the 
chemistry of the H II region and the surrounding neutral shell. 

5. CONCLUSIONS 

We model quasiradial accretion onto seed massive black 
holes in metal and dust poor protogalaxies with the goal of 
evaluating the common assumption that the Mbh ^100 M 
seed black holes accrete at the Bondi accretion rate moderated 
by the Eddington limit. After considering radiative feedback 
effects in the neighborhood of an accreting black hole, we are 
able to derive the following conclusions: 

The local thermal and statistical equilibrium temperature of 
a photoionized gas is a strong function of the metallicity of the 
gas at radii from the black hole where gas becomes captured 
by the black hole. The photoionized gas outside the sonic 
radius is normally in equilibrium, but, particularly for metal- 
licities Z < 0.1 Z Q , the gas passing through the sonic radius 
may not be in equilibrium. For radiative efficiencies e < 0.1, 
the gas will not complete thermal runaway that would heat it 
to the Compton temperature prior to crossing the sonic radius 

Due to the radiative heating of the gas, the sonic radius is 
much smaller than the Bondi radius evaluated in terms of the 
temperature of the protogalactic medium far from the black 
hole; the corresponding gasdynamical accretion rate is thus 
reduced to a small fraction of the Bondi value. The gas tem- 
perature at the sonic radius is also a strong function of the 
radiative efficiency of the accretion near the event horizon of 
the black hole. 

The black hole is surrounded by an H II region, which 
is further surrounded by a photodissociation region if the 
ambient protogalactic medium contains molecules. The ra- 
dius of the H II region in dense (n > 10 5 cm" 3 ) clouds, that 
must be present to support near-Eddington accretion, is small 
(>"ion ~ 1 pc) compared to the size of the protogalaxy. The 
sonic radius is smaller by orders of magnitude than the radius 



of the H II region. This vast separation of scales renders self- 
consistent simulations of accretion onto seed black holes from 
a protogalactic environment particularly challenging, because 
to appropriately capture the radiative effects, the simulations 
must at least resolve the sonic radius. 

Photoionization radiation pressure from the continuum ra- 
diation produced near the event horizon may further diminish 
the accretion rate. Compared to the gravitational accelera- 
tion, the acceleration due to radiation pressure is the strongest 
in the bulk of the H II region, i.e., not too close to the black 
hole where the equilibrium neutral fraction is negligible. We 
have derived a simple criterion for accretion rate suppression 
by the photoionization radiation pressure. 

A significant fraction of the luminosity of the black hole is 
converted into Lya resonance line radiation; resonance line 
scattering keeps this radiation trapped in the neighborhood of 
the H II region. Radiation pressure from the trapped radia- 
tion can exceed the external, confining gas pressure. We do 
not self-consistently solve for a stationary accretion flow in 
the presence of resonance line radiation pressure, but our es- 
timates do suggests that the resonance line pressure could be 
detrimental to stationary accretion at high accretion rates. 

As an alternative to strictly stationary accretion, we out- 
line a model for episodic accretion in which short episodes of 
sudden, rapid accretion alternate with longer periods of accre- 
tion reduced by photoheating and radiation pressure. In this 
model, the outflow and central rarefaction induced by photo- 
heating and radiation pressure acceleration in the bulk of the 
H II region drives down the central accretion rate rapidly. The 
time-average accretion rate is then set by the peak accretion 
rate, divided by the logarithm of the number of e-foldings in 
the outflow-induced accretion rate decay. Since the peak ac- 
cretion rate need not be limited by the thermodynamic and 
kinematic constraints that apply to steady-state solutions, the 
episodic accretion rate may exceed the steady-state rate by a 
large factor. 

Our idealized treatment here already indicates how ex- 
tremely complex the self-consistent accretion problem is. To 
more fully explore the time-dependent accretion flow onto a 
seed black hole, numerical simulations are clearly needed. Of 
particular importance is to study the effect of angular momen- 
tum, the impact of turbulence, and the possible emergence of a 
multi-phase medium in the inf ailing gas. We will report on the 
results from such simulations and on a comparison with the 
analysis presented in this work elsewhere (Milosavljevic et al. 
|2008| Couch et al., in preparation). If our result that the initial 
accretion onto stellar seed black holes is greatly reduced in 
the presence of radiative feedback holds up even in fully three 
dimensional radiation-hydrodynamical simulations, the need 
for massive black hole formation by radiatively-inefficient ac- 
cretion or by direct collapse of massive primordial gas clouds 
might be significantly incre ased (e.g.,|Bromm & Loeb|2003| 
|Begelman et al.|2006|[2008HDjorgovski et al.|2008| >. 
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